Entanglement of formation for symmetric Gaussian states 
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We show that for a fixed amount of entanglement, two-mode squeezed states are those that maxi- 
mize Einstein-Podolsky-Rosen-like correlations. We use this fact to determine the entanglement of 
formation for all symmetric Gaussian states corresponding to two modes. This is the first instance 
in which this measure has been determined for genuine continuous variable systems. 
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One of the main tasks of Quantum Information Theory 
is to quantify the entanglement and the quantum corre- 
lations that quantum states possess. For that, several 
entanglement measures have been introduced during the 
last years Tj. In particular, two of such measures stand 
out for their well defined physical meaning: the entan- 
glement of distillation and of formation (and the cor- 
responding asymptotic generalization, the entanglement 
cost) |2j. They quantify the entanglement of a state in 
terms of the pure state entanglement that can be distilled 
out of it 0] and the one that is needed to prepare it , 
respectively. 

Despite a considerable effort, for the moment we can 
only evaluate the entanglement of formation (EoF), or 
the entanglement of distillation for a few sets of mixed 
states. The reason is that these quantities are denned 
in terms of an optimization problem which is extremely 
difficult to handle analytically. Despite this fact, in a 
remarkable work Wootters managed to derive an an- 
alytical expression for the EoF for all two-qubit states. 
The EoF has been also determined for highly symmetric 
states (isotropic states p, and Werner states jg). These 
expressions are important theoretical tools. From a more 
practical point of view, they can be applied to quantify 
the entanglement created in current experiments as well 
as to compare the capability of different experimental 
set-ups. For low dimensional systems without symme- 
tries, one can still use numerical methods to determine 
the EoF although they are often not very efficient. 
For infinite-dimensional systems, however, a numerical 
approach is not feasible. 

Among all quantum states in infinite dimensional sys- 
tems, Gaussian states play an important role in quantum 
information. From the experimental point of view, they 
can be created relatively easily || , and one can use them 
for quantum cryptography [lfj and quantum teleporta- 
tion 0, 0|. On the theoretical side, separability 
and distillability 0] criteria for bipartite systems have 
been fully developed. Moreover, pure Gaussian states 
are intimately related to Heisenberg's uncertainty rela- 
tion since they minimize such a relation for position and 
momentum operators. 



In this work we determine the EoF of all symmetric 
Gaussian states of two modes. Those states arise nat- 
urally in several experimental contexts. For example, 
when the two output beams of a parametric down con- 
verter are sent through optical fibers or in atomic 
ensembles interacting with light [l^. In order to de- 
termine the EoF, we connect the entanglement of pure 
states, as measured by the von Neumann entropy of the 
restriction, with the kind of correlations established by 
Einstein, Podolsky and Rosen (EPR) in their seminal pa- 
per 0. In fact, we show that two-mode squeezed states 
|l7| play a very special role in this relation, since they are 
the least entangled states for a given correlation of this 
type. This provides a new characterization of two-mode 
squeezed states. Finally, we show that the decomposi- 
tion that leads to the EoF is a decomposition in terms of 
Gaussian states. 

We consider two modes, A and B, with corresponding 
Hilbert space H — Ha ® Hb and canonical operators 
Xa,b and Pa,b- The two-mode squeezed states have 
the form 
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where r > is the squeezing parameter, and \N) denotes 
the N-th Fock state, i.e. a^a\N) A = N\N) A , tfb\N) B = 
N\N) B , where a — (Xa + iPA)/V2 and b = {Xb + 
zPb)/ 1 /^ are annihilation operators. 

In the following, we will denote by ip an arbitrary nor- 
malized state in H. We define its EPR-uncertainty as 
follows: 



A(V0 



X T 



Al(P A + P B ) 



,(2) 



where, as usual 



(Xip\X>ip) - (V>|AT|V>) 2 , setting 
A 2 (A) = oo if ip is not in the domain of A. Clearly, 
A(V>) £ (0,1]- This quantity measures the degree of 
non-local correlations, and would be zero for the ide- 
alized state considered by Einstein, Podolsky, and Rosen 
For any state, A(^) < 1 implies the existence 
of such non-local correlations. Note that this condi- 
tion is met only if at least one of the uncertainties of 
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(X A - X B )/V2 or (P A +P B )/V2 lies below 1 (the stan- 
dard quantum limit). This implies that the correspond- 
ing states must possess a certain squeezing. In fact, the 
two-mode squeezed states are standard examples of 
states displaying these correlations since 



A[*,(r)] = e- 2r < 1. 



(3) 



Any value of A 6 (0, 1) is achieved by the two-mode 
squeezed state with squeezing parameter 



r A :=--ln(A) 



(4) 



The EPR-uncertainty of a given state tp is certainly 
related to its entanglement. For pure states this last 
property is uniquely quantified by the entropy of entan- 
glement, E(ip), which can be determined as follows. Let 
us write the Schmidt decomposition of ip as 



^ c n \u n )a ® \vn)b, 



(5) 



N=0 



{itjv} and {vjy} are orthonormal bases in H.a,b, respec- 
tively and c — (c , Ci, . . .) G C where 



C:={cell 



l,cjv>c*r+i>0VJV}. 



Then 



EM = e(c) 



£ 4iog(4). 

jv=o 



(6) 



Note that this quantity can be infinite for some states. 
For the two-mode squeezed states (QJ we have 

£[f s (r)] = cosh 2 (r) log[cosh 2 (r)] -sinh 2 (r) log[sinh 2 (r)]. 

With the above definitions we can state the special 
role that two-mode squeezed states play in relation with 
EPR-correlations and entanglement: 

Proposition 1: For all tp e H, E(ip) > E[^ s (r A ^)}. 

In order to give a clear interpretation of this result, 
we reformulate it in two equivalent forms: (i) For any 
given A G (0,1), E[9 a (r A )] = inf^{^(V)} with V 
fulfilling A(V>) = A; (ii) For any given E G (0,co), 
A[* a (r)] = inf 1> {A(V')} with V fulfilling E(ip) = E 
and r such that E[^ s (r)] — E. Note that the equiv- 
alence of this last formulation is ensured by the fact 
that £ , [^I / S (r)] and A[^ s (r)] are monotonically increas- 
ing and decreasing functions of r, respectively. The 
first statement characterizes two-mode squeezed states 
as the cheapest (regarding entanglement) to achieve a 
prescribed EPR uncertainty. The second one character- 
izes them as those with maximal EPR-correlations (min- 
imal A) for any given value of the entanglement. 

In order to prove Proposition 1 we introduce two lem- 
mas and the following definition. Given ceCwe define 



6(c) := 1 + 2 (4 - c N c N -x)N. 



(7) 



We have 6(c) < 1 and 6(c) = A(ip) whenever \un) — 
\vn) = | AO [cf. ©]. 

Lemma 1: For all ip with Schmidt decomposition JSJ, 
A(ip) > 6(c). 

Proof: Since 6(c) < 1 we can restrict ourselves to ip 
with A(tp) < 1. Without loss of generality we can assume 
that (ip\a\x/)) — (tp\b\tl>) = 0. Otherwise we can always 
find ip 1 fulfilling this condition, with the same Schmidt 
coefficients as i/j and with A(ip') — A(ip) [2oj |. We have 

oo 

a^a\uN) + (Vjy\tfb\Vjy)) 

N=0 
oo 

- Y c N CM((uN\a\u M )(vN\b\v M ) + c.c). 

N,M=0 

> xnin[Z(u), Z(v)], 

where 



Z(u) := 1 + 2 c 2 N (u N \a^a\u N ) 

oo 

-2 Y c N c M \(u N \a)\u M )\ 2 - 

N,M=0 

Without loss of generality let us assume that 
mm[Z(u), Z(v)] = Z(u) =: Z. We can rewrite it as 
Z = En=oSm=w+i( Cw — cm) Xn^m, where Xjs[ y M '■= 
|(ujv|a^|itM)| 2 + KumMIujv)! 2 - Now, since c G C we can 
write (cat - c M ) 2 > J2r=n( c R ~ c R+i) 2 , for M > N + 1, 
so that 



Z > 5>«-c« +1 ) 2 ]T J2 X " 

R=0 N=0 M=R+1 

oo 

= ^( C fl- C7?+1 ) 2 (i? + i + 2r fl ), 



where 



R=0 



{u N \a' N a N \u N ) 
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with o/v := a— (ujsr\a\uN) . Now, using that un -L um for 
N ^ M we have (ujv|ci^|um) — (un\o>n\um) which, to- 
gether with J2a=M^N \(un\^ n \u m )\ 2 < (u N \a^ N a N \u N ) , 
yields that Yr > for all i? and therefore 



oo 

Aty) > Z > E ( c * " cfl+i) 2 (i? + 1) > 6(c). 

R=0 



□ 



JV=0 



Lemma 1 indicates that for a given set of Schmidt co- 
efficients c E C EPR-correlations are maximized if the 
Schmidt vectors are chosen to be Fock states in the right 
order, i.e. |itjy) = \vn) — \N). Next we will show that 
for fixed A, the choice of Schmidt coefficients minimizing 
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the entropy of entanglement is given by those of a two- 
mode squeezed state. Since the entropy and the EPR 
entanglement are explicitly known functionals e(c) and 
6(c) on the sequences c £ C, this is a classical constrained 
variational problem. 

Lemma 2: For A £ (0, 1), and any sequence c € C 
with <5(c) = A, we have e(c) > e(c A ) = E[^ s (r A )], 
where c N oc exp(— Nr A ) is the unique geometric sequence 
in C with <5(c A ) = A. 

Sketch of proof: We apply the method of Lagrange 
multipliers for constrained minima to the infinitely many 
variables cq, c\ . . ., leaving aside the technicalities of mak- 
ing this rigorous. These involve restricting c to finite di- 
mensional spaces, then letting the dimension of the space 
tend to infinity, and controlling the attained minima in 
this limit. 

With a choice of Langrange multipliers fi and A > 0, 
designed to simplify the expressions to come, we are thus 
looking for stationary values c £ C of the functional 

F(c,X,,) : =e(c) + ^[«(c)-A] + ^(||c||-l). 
We obtain 

2c N [NX + n _ = \[N c/v— i + {N + l)c N+1 }. (8) 

where we have defined c_j = 1. One can immediately 
see that cn > and thus we can divide JHJ) by cn and 
subtract the same expression but for N + 1. Defining 



x N := c N+ i/c N 



-2r N 



£ (0,1] for N = 0,1,. 



writing A = 2r / sinh (r) for some r > we find 
£at+i = x n — A N — B N , 



and 



(9) 



where N = 0, 1, . . . and 



B N 



N + 2 

N 
N + 2 



sinh 2 (r N ) - —sinh 2 (r) , (10a) 
r 



1 



1 



X N XJV-1 



(10b) 



If we fix r > and xo, we have three possibilities: (i) 
xo < e~ 2r . Then, by induction, xjv is decreasing, and will 
reach some xjv < for finite N, which is impossible; (ii) 
xq > e~ 2r . Then xn is increasing, and the normalization 
condition for c cannot be fulfilled. Hence we must have 
the third possibility (iii) xq — e~ 2r , which implies that 
xjv = e~ 2r for all N. Hence cat is a geometric sequence 
oc exp(— 2iVr).m 

With the help of Lemmas 1 and 2 we are now in the 
position of proving Proposition 1. 

Proof of Proposition 1: Given i\> £ H, if A(W) — 1 then 
it is trivial. Otherwise, using (j^J and Lemma 2 we have 

Ety) = e(c) > E[* s {r s{c) )\ > E[V s (r AW )}, (11) 

where for the last inequality we have used ?"a(j/>) ^ r S(c) 
[which follows from Lemma 1 and ©] and the fact that 
E[$! s (r)] increases monotonically with r. p | 



In the following, we will apply Proposition 1 to deter- 
mine the EoF of symmetric Gaussian states of two modes. 
For a given density operator, a, we define its covariance 
matrix (CM) 7 as usual, 

7y - := ix[(IURj + RjRi)p] - 2tT(R i p)tr(R j p), (12) 

where {R t ,i = 1,..,4} := {X A , P A , X B , P B }. Up to lo- 
cal unitary operations, it can always be written in the 
standard form 



/ n k 2 

n 

k x m 

\ — A: p m 



\ 

■fe 



(13) 



We will concentrate here on symmetric states, i.e. those 
which are invariant under exchange of subindices A and 
B and therefore fulfilling m = n. Without loss of gen- 
erality we can choose k x > k p > 0. In this case, 7 is a 
CM iff n 2 — fc 2 > 1 and describes an entangled state iff 
1 > (n-k x )(n-k p ) E3. Next we apply local (unitary) 
squeezing transformations to yield the state in a more ap- 
propriate form without changing its entanglement prop- 
erties. In the Heisenberg picture, the transformation 
multiplies (divides) Xa,b (Pa,b) by [(n~k p )/(n~k x )] 1 ^ 4: . 
A simple calculation gives 



A(cr) = \J(n- k x )(n - kp) =: S 



(14) 



where A(er) is defined analogously as in |2"|l. 

Our goal is to determine the EoF of a. This is defined 
as Ep(a) := inf £>£(£>), where the infimum is taken with 
respect to all sets of the form D = {pk,ipk} which give 
rise to a decomposition of cr, i.e., 



o- = ^2pk\*l>k)(ipk 



(15) 



where the tpk £ 7~L are normalized and pk > 0. Note that 
the sum can run over continuous indices. For the set D 
we define 



£(D) :=J2PkE^k)- 



(16) 



We call the set D a decomposition of a. A particular 
decomposition Dq of a is defined through 



a oc 



d£W($\9.(rs)){9.(r s )\W0e 



where W(£) = e 1 ^ R is the Weyl displacement operator 
and 75 < 7 is the CM of the two-mode squeezed state 
with squeezing parameter r<j. Since W^(x) are local 
unitary operators, we have £(D n ) = E[^ s (rs)]. 

We also introduce the auxiliary function / : (0, 1] — » 
[0,oo) 



/(A) 



_(A)log[c + (A)]- c _(A)log[c_(A)], (17) 



4 



where c±(A) := (A -1 / 2 ± A 1 / 2 ) 2 ^. One can readily 
show that / is a convex and decreasing function of A 
and that 



E[* s (r A )] = f(A). 



(18) 



Proposition 2: Ep(a) = f[-J (n — k x )(n — k p 



Proof: We just have to prove that for any decompo- 
sition D, £{D) > f(5), where 5 is given in (|14|l . since 
the decomposition Do already achieves this value, i.e. 
£{D ) = E[9 a (r s )] = f(S) [c.f. Eq. For any de- 

composition we have 



£(D) > 5>*/[A(Vfc)] > f 



The first inequality is a consequence of Proposition 1 
and l|18|l . The second is due to the convexity of /. 
Finally, the last one is a consequence of the fact that 
5 > ^pfeA(V'fc) (which can be easily checked by using 
the Cauchy-Schwarz inequality) together with the fact 
that / is a decreasing function of its argument, rj 

In summary, we have determined the EoF of symmet- 
ric Gaussian states by establishing a connection between 
EPR-like correlations and the entanglement of a state. 
The result implies that the measured quantities in some 
of the recent experiments dealing with atoms 0] and 
photons |2l| not only qualify entanglement but also quan- 
tify it. We expect that the methods introduced here 
will allow to determine the EoF and other properties of 
more general Gaussian states. The optimal decomposi- 
tion Do that gives rise to the EoF is a mixture of Gaus- 
sian pure states, which means that those states are the 
cheapest ones in terms of entanglement to produce sym- 
metric Gaussian states. Thus, it is tempting to conjec- 
ture that this is also true for all Gaussian states. Finally, 
the results presented here provide a new characterization 
of two-mode squeezed states as those, which achieve a 
maximal EPR-like correlation for a given value of the 
present entanglement. 
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